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Abstract: We introduce a new model for the dynamics of centroblasts and centrocytes in a 
germinal center. The model reduces the germinal center reaction to the elements considered 
as essential and embeds proliferation of centroblasts, point mutations of the corresponding 
antibody types represented in a shape space, differentiation to centrocytes, selection with 
respect to initial antigens, differentiation of positively selected centrocytes to plasma or 
memory cells and recycling of centrocytes to centroblasts. We use exclusively parameters 
with a direct biological interpretation such that, once determined by experimental data, 
the model gains predictive power. Based on the experiment of Han et al. (1995b) we predict 
that a high rate of recycling of centrocytes to centroblasts is necessary for the germinal 
center reaction to work reliably. Furthermore, we find a delayed start of the production 
of plasma and memory cells with respect to the start of point mutations, which turns out 
to be necessary for the optimization process during the germinal center reaction. The 
dependence of the germinal center reaction on the recycling probability is analyzed. 
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1 Introduction 



Germinal centers (GCs) are essential elements of the humoral immune response (for a 
review see MacLennan, 1994). In a primary immunization some B-cells respond to the 
presented antigen and are activated. These active cells migrate to follicular dendritic 
cells (FDCs) present in lymphoid tissues. Unprocessed fragments of the antigen, so called 
epitopes, are expressed on the FDCs. In the milieu of the FDCs the B-cells are in a phase of 
intense proliferation and show the essential features of B-blasts such that they are denoted 
as centroblasts. These proliferating centroblasts together with the FDCs develop to a GC. 

This morphological development is initiated through an unidentified differentiation sig- 
nal after about three days. The resulting GC is characterized by two zones (Nossal, 1991): 
the dark zone filled with continuously proliferating and mutating centroblasts and the light 
zone containing the FDCs and centrocytes. The latter are generated from the centroblasts 
(Liu et al., 1991), which are available in great diversity (Weigert et al., 1970) because of 
a very high mutation rate (Nossal, 1991; Berek & Milstein, 1987). In contrast to cen- 
troblasts the centrocytes express antibodies of cell specific types on their surface. In this 
way they are able to undergo a specific selection process through the interaction with the 
epitopes presented on the FDCs in the light zone. The diversity of the proliferating and 
mutating centroblasts together with the specificity of the selection process in the light zone 
allow an optimization of the affinity between the antibodies and the epitopes. This affinity 
maturation in the GC leads to a dominant population of cells with high affinity of the 
corresponding antibodies to the epitopes within 7-8 days after immunization (Jacob et al., 
1993). 

Centrocytes undergo apoptosis if they are not positively selected in an appropriate 
range of time. They are rescued from apoptosis (Liu et al., 1989; Lindhout et al., 1993; 
Lindhout et al., 1995; Choe & Choi, 1998) if they are able to bind to the antigen on the 
FDCs and to interact successfully with T-helper cells present in the external region of the 
light zone. Even antigen independent signals may be required for the centrocytes to survive 
(Fischer et al., 1998). Successfully selected centrocytes receive a differentiation signal which 
determines their fate. They may differentiate into memory cells that are important in the 
case of a second immunization. Alternatively, they differentiate to antibody producing 
plasma cells. These are intensifying the immune response in two ways: the number of 
antibodies increases and their specificity with respect to the direct unspecific immune 
response without GC reaction is optimized. In both cases the corresponding cell is leaving 
the GC environment. 

A third possibility exists for the positively selected centrocytes: they may be recycled 
to centroblasts and reenter the highly proliferating and mutating stage of development 
in the dark zone (Kepler & Perelson, 1993). This recycling hypothesis has been neither 
directly checked experimentally (for an indirect check see Han et al., 1995b) nor - if true 
- the probability of recycling is known. It would be of great interest to know about the 
necessity of a centrocyte recycling process for the affinity maturation. The most evident 
argument to support the recycling hypothesis is based on the occurrence of multi-step 
somatic hypermutations in the GC reaction. Under the assumption of random hypermu- 
tations, recycling of already mutated cells to fast proliferating centroblasts appears to be 
a necessary process. 
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In this paper we estimate how large the recycling probability should be, i.e. which 
proportion of the positively selected centrocytes should be recycled to get a reliable GC 
reaction. For this purpose we introduce a simplified model of the GC dynamics based on 
a functional view on it. The model differs from others (Oprea & Perelson, 1996; Oprea & 
Perelson, 1997; Rundell et al., 1998) by omitting some details of the GC reaction which we 
consider as unnecessary to understand the obligatory occurrence of the recycling procedure. 
On the other hand compared to the model introduced in (Kesmir & de Boer, 1999) we 
added a formalism to treat somatic hypermutation in a shape space, which seems to be 
essential to understand the process of affinity maturation. A very important difference to 
related models is that we embed the experiment (Han et al., 1995b) for the first time. This 
experiment provides an indirect evidence for the existence of the recycling process and we 
will interpret it in a quantitative way in order to extract information on the recycling 
probability from it. We use exclusively parameters with a direct biological interpretation 
to ensure the predictive power of the model. 

In the next section we will define the shape space, in which the cell population dy- 
namics takes place, and the corresponding dynamical quantities and parameters including 
their dynamical interdependence. All parameters introduced in the model are determined 
by experimental observations. The third section is dedicated to the presentation and in- 
terpretation of the model in comparison with experimental GC reactions. We will analyze 
the stability of our results to verify the notion of a prediction. Finally, we investigate the 
dependence of some characteristics of the germinal center reaction on the variation of the 
recycling probability and give some concluding remarks. 

2 Definition of the model 

2.1 Shape space and affinity space 

Postulates Before defining a model for the GC dynamics it is necessary to specify the 
space in which the dynamics takes place. As we want to represent cells corresponding to 
different types of antibodies we use the well known shape space concept (Perelson & Oster, 
1979). This ansatz is based on the following assumptions: 

• Antibodies can be represented in a phenotypical shape space. 

• Mutation can be represented in the shape space by next-neighbor jumps. 

• Complementarity of antibodies and antigen. 

• Homogeneity of the affinity weight function on the shape space. 

Representation of antibody phenotype The shape space is taken as a D-dimensional 
finite size lattice of discrete equidistant points, each of them representing one specific 
antibody type. Thus the shape space becomes a phenotype space, i.e. it is not primarily a 
representation of genetic codes but of the resulting features due to specific genetic codes. 
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Representation of mutation Nevertheless, we want to represent point mutations de- 
fined in an unknown genotype space in the shape space. To this end we are compelled to 
formulate the action of a point mutation in the shape space: 

A point mutation of a given antibody is represented in the phenotypical shape space by a 
jump to one of the nearest neighbor points. 

This assumption is not very spectacular as it only requires that a point mutation does 
not lead to a random jump of the antibody phenotype, i.e. that conformation, electrical 
properties, etc. are not dramatically changed. Surely, the mutation of some key base pairs 
may exist which imply fundamental changes of the encoded antibody features. However, 
these exceptions will not alter the dynamics of the GC until the number of such key base 
pairs remains small with respect to the number of smooth mutation points. 

Antigen representation and complementarity A central quantity is the affinity of a 
given antibody and an antigen epitope. Having a representation of antibodies in the shape 
space, a counterpart for antigens is necessary. We emphasize that the number of possible 
antibodies is finite, whereas the diversity of antigens is principally unbounded. Therefore, 
we represent an antigen on the B-cell shape space at the position of the B-cell with highest 
affinity to the antigen, i.e. at the position of an antibody with optimal complementarity to 
the antigen. 

Property of smooth affinity Starting from these assumptions (phenotype shape space, 
representation of mutation by next-neighbor jumps, and complementarity) we have reached 
a representation of antibodies and antigens in a shape space, which has the property of 
affinity neighborhood, i.e. that neighbor points in the shape space have comparable affinity 
to a given antigen. This property is a direct consequence of the definition of mutations 
in a (phenotype) shape space. Starting from a non-directed mutation as base for affinity 
maturation in GCs, we estimate the affinity neighborhood as a necessary property of the 
underlying shape space. An initial B-cell with a particular affinity to a certain antigen 
must have the possibility of successively optimizing the affinity, i.e. of stepwise climbing 
an affinity hill. If no affinity neighborhood would exist, i.e. if there was a random affinity 
distribution on the shape space, each mutation would lead to a random jump in the affinity. 
An optimization of the affinity to the antigen may occur accidentally in this scenario, but 
should be a very rare event. 

Homogeneous affinity weight function Affinity neighborhood allows for the defini- 
tion of an affinity weight function which determines the affinity of antigen and antibody in 
dependence of their distance in the shape space. We assume this weight function a(0, 0*) 
to apply equally well in all regions of the shape space which corresponds to a homogeneous 
affinity distribution over the shape space, and to be of exponential form: 



where V is the width of the affinity weight function and 1 1 • 1 1 denotes the Euclidean metric. 
The question remains with which exponent i] the distance enters the exponential weight 




(1) 
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function and it will be argued that 77 = 2 is a reasonable value (see Fig.|2|). 
2.2 Formulation of the dynamics 

Having a shape space at hand it is possible to define distributions -8(0), C(<j>), A(4>), and 
0((f>) of centroblasts, of centrocytes, of the presented antigen epitopes, and of the plasma 
and memory cells on the shape space (see Tab. [I]). We focus on the centroblast distribution 
B and analyze its dynamical behavior in the different phases of the GC reaction. 



Cell type 


Variable 


Processes 


antigen 


m 


interaction with centrocytes 


centroblasts 


5(0) 


proliferation 






mutation 






differentiation to centrocytes 


centrocytes 


c{<t>) 


selection and apoptosis 






recycling to centroblasts 






differentiation to plasma and memory cells 


plasma and memory cells 


0(0) 


removed from GC 



Table 1: The variables in the mathematical model of the GC reaction and the processes that 
are considered. 



There were several attempts to divide the germinal center reaction in different working 
phases. From the point of view of our model we are led to a new functional phase distinction 
of the GC reaction (see Fig. |I]). Our three-phase description of the GC reaction is in 
accordance with most of the models established so far (see e.g. MacLennan, 1994; Liu et 
al., 1991; Kelsoe, 1996) but is in contradiction to the time scales found in (Camacho et a!., 
1998), which differ strongly. It is not intended to find mechanisms of transition between 
the phases of the germinal center reaction. These dynamical phases are assumed according 
to experimental observations. 

The first phase is the already mentioned phase of fast proliferation of a low number 
of seeder cells in the environment of FDCs. In this phase it is likely that somatic hyper- 
mutation is not taking place to a relevant amount (Han et al., 1995b; Jacob et al., 1993; 
McHeyzer- Williams et al., 1993; Pascual et al., 1994b) such that it may be understood as 
an enlargement phase of the cell pool available for the following process of affinity matu- 
ration. The corresponding dynamical behavior of the centroblast distribution is described 
as 

dB 

— (0) = pB{(j>) for t-t <Ah , (2) 

where p denotes the proliferation rate and to is the time of first immunization. This phase 
lasts At 1 = 3 days and is well established by experiment (Liu et al., 1991). 

After three days the optimization phase starts. The GC gets its morphological form, 
i.e. dark and light zones emerge. In the dark zone centroblasts continue to proliferate but, 
additionally, somatic hypermutations are broadening the initial centroblast distribution on 
the shape space. At the same time the selection process operates in the light zone and 
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Proliferation (p) 

Mutation (m) 
Selection (<2q) 

Differentiation (g) 
Recycling (1-q) -1? 

Plasma and memory 
cell production (q) — »~ 



t Q fj t 2 t 3 

At, ► ► 

i n in 

Figure 1: The three model phases of the germinal center reaction: (I) proliferation phase 
(starting at t = —3d), (II) optimization phase (starting at t\ = Od, the duration A.t 2 is 
determined in Sec. |2.8| ), (III) output phase (ending at t 3 = 18 d). The parameters in brackets 
refer to the model equations Eq. (|2|), Eq. (^]), Eq. (^|), and Eq. (|8|) and can be determined from 
experimental data (see Sec. ^2] - [2.8|) . 

the selected centrocytes are either recycled to centroblasts or differentiate to plasma and 
memory cells. The development of the centroblast distribution on the shape space is now 
described by 

^(0) = (p- 2pm -g)B(<j>) + ^ £ Z?(0 + A</>) 

dt D ||A0||=1 

+ (l-g)^C7(0)A(0*)aoa(0,0*) , (3) 

for t — to > Ati. Here, g is the rate of centroblast differentiation into centrocytes, m is 
the mutation probability, a is the probability of an optimal centrocyte to be positively 
selected, q is the probability of differentiation into plasma and memory cells for positively 
selected centrocytes, and D is the dimension of the underlying shape space. The different 
contributions to Eq. @ are discussed in detail below. Note that this equation does not 
include eventual finite size effects due to small populations. 

It is not clear a priori, if differentiation into antibody producing plasma and memory 
cells is triggered already in this second phase. To allow a start of the output cell production 
delayed by the time interval At 2 we divide the optimization phase into two sub-phases 
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which differ in the output probability q: 



q = for < t - At x < At 2 
q > for < t - At 2 < At 3 



(4) 




The time delay At 2 will be fixed by experimental data (see Sec. |2.8|) . The output phase in- 

clude optimization and production of plasma and memory cells and lasts for the remaining 
GC life-time, which is about At 3 = 21days-3days- At 2 (Choe & Choi, 1998; Kelsoe, 1996; 
Jacob et al., 1991). Then, only a few proliferating B-blasts remain in the environment of 
the FDCs (Liu et al., 1991). 

Proliferation 

During the whole GC reaction a fast proliferation of B-cells takes place. After the activation 
of B-cells by an interaction with an antigen, these move to FDCs and undergo a fast 
proliferation phase in their environment (Liu et al., 1991; Hanna, 1964; Zhang et al., 1988; 
Grouard et al., 1995). It is likely that the fast proliferation is triggered by the dendritic 
cells (Dubois et al., 1999). After about three days (Liu et al., 1991) when the GC starts to 
develop its characteristic light and dark zone (Camacho et al., 1998), a fast proliferation 
of centroblasts is still observed in the dark zone. The extremely high rate of proliferation 
could be determined to be (Liu et al., 1991; Zhang et al., 1988) 



a result known already for a long time (Hanna, 1964). The population of centroblasts at 
grows exponentially, which is represented by the term pB((p) in Eq. (0) and Eq. ([3]). 



It is likely that somatic hypermutation does not occur in the proliferation phase of the GC 
reaction (Han et al., 1995b; Jacob et al., 1993; McHeyzer- Williams et al., 1993; Pascual 
et al., 1994b). On the other hand the growth of the centroblast population is reduced by 
possible mutations from state to its neighbors + A0, where | |A0| | = 1. In a continuous 
space this corresponds to a diffusion process as used in (Perelson & Wiegel, 1999) to 
represent mutation. Each pair of cells produced by proliferation will populate a neighbor 
state with a probability of m. This mutation probability turns out to take extremely high 
values of m ~ 1/2 (Nossal, 1991; Berek & Milstein, 1987), which corresponds to a factor 
10 7 larger probability compared to mutations outside the GCs (Janeway & Travers, 1997). 
We want to point out that here only point mutations of phenotypical relevance are taken 
into account. As a consequence the population growth by proliferation is reduced by the 
important amount 2pm, which results in an effective proliferation rate at of p(l — 2m). 
One observes that for m = 1/2 the two first terms in Eq. (||) cancel. This corresponds to 
the situation that in each centroblast replication one new cell remains in the old state and 
the second new cell mutates to a neighbor state such that the total number of cells in state 
remains unaltered. 



Mutation 
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In the same way as mutation depopulates the state in the shape space and populates 
its neighbors, is populated by its neighbors. This gives rise to the non-diagonal term in 
Eq. ([!]). Each neighbor of mutates with the same rate 2pm from + A0 to one of its 
2D neighbors. Only the mutation from + A0 to the particular neighbor enhances the 
population of centroblasts at such that the rate is weighted by the inverse number of 
neighbors 1/(2D) and is summed over all possible neighbors of 0. 

Selection, recycling and cell production 

After the establishment of the light and dark zones in the GC the differentiation of cen- 
troblasts to antibody presenting centrocytes is triggered to allow a selection process in the 
light zone. The centroblasts in the shape space state are diminished with the rate of 
differentiation g, leading to the term — g£>(0) in Eq. (||). The centrocyte population C(0) 
of type is enhanced simultaneously by the same amount of cells: 

C(0) = +gB(<P) . (6) 

The centrocytes move to the light zone where they undergo a selection process. Their 
further development is splitted three-fold. The non-selected centrocytes die through apop- 
tosis and are eliminated from the GC dynamics. It is known for a long time that the 
centrocytes undergoing apoptosis were in cell cycle a few hours ago (Fliedner, 1967) such 
that it is likely that they differentiated from the centroblasts. On the other hand apoptosis 
takes place in the light zone (Hardie et al., 1993) where centroblasts are not present in 
high density. The selected centrocytes are emitted from the GC with probability q and dif- 
ferentiate either into antibody producing plasma cells or to memory cells. The model does 
not distinguish between plasma and memory cells. Only their sum is taken into account 
and is denoted by output cells 0(0). Nevertheless, the dynamics of both types of output 
cells may be different (Choe & Choi, 1998). Also, the degree of affinity maturation differs 
(Smith et al., 1997). In addition, we do not consider any further proliferation of output 
cells in or outside of the GC, which may be important for a quantitative comparison with 
experimental measurements. 

Alternatively the selected centrocytes are recycled to become centroblasts and to reen- 
ter the proliferation process in the dark zone. This happens with probability 1 — q and 
contributes to the centroblast distribution, corresponding to the last term in Eq. ([|). 

We want to emphasize that we do not resolve the time course of the selection process, 
which is regarded to be fast with respect to the centroblast proliferation time scale. This 
procedure is justified by a minimal duration of a typical selection process of 1 — 2 hours 
(van Eijk & de Groot, 1999), which is about one fourth of the proliferation time scale. 
Nevertheless, one should be aware of possible implications due to the fact that the selection 
process does not occur instantaneously. We effectively enclose the selection time in the 
parameter g, which in this way becomes a measure of a complete selection cycle including 
the differentiation into centrocytes, a first selection at the FDCs, a second selection at 
T-helper cells and finally the recycling process. This has two consequences for the model: 
The number of selectable centrocytes for each shape space state is given at every time 
by Eq. (^), i.e. the number of centrocytes just being in the selection process. Furthermore, 
the details of the multi-step selection process (Lindhout et al., 1997) are omitted. The 
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selection is modeled by a sum over the shape space of the antigen distribution presented 
on the FDCs weighted by the affinity function Eq. (fj). In other words the centrocytes 
gB(<f)) are selected if an antigen is close enough in the shape space. The meaning of close 
enough is determined by the width of the affinity function and by the probability Oq of a 
positive selection for an optimal centrocyte with respect to one presented antigen. 

One should be aware that the probability oq is not necessarily equal to one, as centro- 
cytes are in the state of activated apoptosis (Cohen et al., 1992) such that their life time is 
finite and they have to be selected within this life time to be rescued from apoptosis. This 
maximum selection probability is determined by the relation of the centrocyte life time of 
about 10 h (Liu et al., 1989, 1994) and the duration of the selection process of 1 — 2 h (van 
Eijk & de Groot, 1999). A rough estimate using two Gaussian distributions with reason- 
able widths for the life time and the selection time peaked at 10 h and 1.5 h respectively 
leads to Oq ~ 0.95. Thus, if a GC contains only centrocytes of optimal complementarity 
to a presented antigen, still about 5% of them will undergo apoptosis. 

In conclusion, selection can be described by the convolution term 

£C(0)A(0> o a(0,0*) , (7) 

<t>* 

where the number of selectable centrocytes is given by Eq. ([]) and A(<p*) represents the dis- 
tribution of antigens on the shape space. This selection term contributes with probability 
q to the production of plasma and memory cells 

^(0) = g£CX0)A(0> o a(0,0*) . (8) 

All selected centrocytes not emitted from the GC are recycled and in this way enhance the 
centroblast population giving rise to the last term in Eq. (|3|). 

2.3 Initial conditions 

Antigen distribution 

Throughout the paper the antigen distribution is assumed unequal to zero at exactly one 
point in the shape space: 

A(<f>*) = A(fa) = 5(0* -<M , (9) 

where S(-) is one for its argument equal to zero and zero otherwise. Consequently, the sum 
in the convolution term in Eq. (Q) is reduced to one contribution. This implies that we 
consider only one type of antigen epitope to be present in the GC. What seems to be a 
restriction at first sight turns out to be an assumption without any loss of generality. In 
view of about 10 11 possible states in the naive repertoire of humans (Berek & Milstein, 
1988) two randomly chosen antigen epitopes will not be in direct neighborhood of each 
others in the shape space. Both antigen epitopes will only have an influence on each other 
during the GC reaction when centroblasts corresponding to one of the antigen epitopes have 
a non-negligible affinity to the other epitope. In the shape space language this situation 
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corresponds to an overlap of two spheres centered at both antigen epitopes. The radius of 
these spheres is determined by the maximum number of mutations which may occur during 
a GC reaction. This number will rarely exceed 9 (Kuppers et al., 1993; Wedemayer et al., 
1997) such that in a 4-dimensional shape space (a more detailed discussion of the shape 
space dimension follows in Sec. |2.5| ) one sphere of this size covers about 10 _8 th part of 
the space. An overlap is unlikely as an unreasonable large number of 10 8 antigen epitopes 
is needed to get an important probability for a mutual influence. Therefore we consider 
one-antigen epitope GC reactions only. A multiple-antigen or a multiple-antigen epitope 
GC reaction has to be considered as sum of one-antigen epitope GCs. 

There exists experimental evidence that the amount of presented antigens is reduced 
during a GC reaction (Tew & Mandel, 1979). Since the amount of antigen is only halved 
during 30 days (Tew & Mandel, 1979; Oprea et al., 2000) we do not expect an important 
recoil effect on the selection probability, which would become time dependent otherwise. 
Furthermore, we believe that this antigen diminution is too weak to be responsible for the 
vanishing cell population in GCs after about 3-4 weeks. Anyhow, we aim to show that 
our three-phase model allows for a vanishing population solely due to the dynamics of the 
system, without the inclusion of any antigen distribution dynamics. 

B-cell distribution 

The B-cells which have moved to FDCs are in a fast proliferation phase. In this first phase 
of GC development the cell population grows exponentially until it has reached about 10 4 
cells after three days. One may ask how many seeder cells are necessary and sufficient 
to give rise to this large cell population within such a short space of time. One finds 
experimentally that the follicular response is of oligoclonal character (Liu et al., 1991; 
Jacob et al., 1991; Kuppers et al., 1993; Kroese et al., 1987) and the number of seeder cells 
is estimated between two and six. This result is in accordance with the proliferation rate 
of ln(2) / (6 h) cited above. 

The number of seeder cells J2^B((p,t = 0) determines the initial B-cell distribution 
B((f)) in our model. Corresponding to the experimental results cited above, we start with 
the reasonable number of three seeder cells for the follicular response to an immunization 
with a specific antigen. These cells are chosen randomly in the shape space within a sphere 
of radius Rq around the presented antigen. Due to the assumed affinity neighborhood this 
reflects a threshold affinity of antigen and activated initial B-cells. This was already found 
experimentally (Agarwal et al., 1998): The activation of a B-cell by an antigen is necessarily 
connected with a minimum affinity of the corresponding epitopes and antibodies. The 
radius of the sphere can be determined by the maximum number of mutations occurring 
during the process of affinity maturation in GCs: 

#o = AUx , (10) 

which is found to be 8 or 9 (Kuppers et al., 1993; Wedemayer et al., 1997). In this picture, 
B-cells with an affinity above threshold are activated by a specific antigen, giving rise to 
a fast immune response with antibodies of low but non-vanishing affinity. The germinal 
center reaction has to be understood as an optimization process for these initially activated 
B-cells leading to a second slower primary immune response of high specificity. 
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2.4 Width of affinity function 

Due to affinity neighborhood in the shape space the affinity function may be chosen ac- 
cording to Eq. (H). The width T has to be consistent with the affinity enhancement that 
is achieved during a GC optimization process. Let the factor of affinity enhancement be 
1CK and the number of corresponding somatic hypermutations be N. Then from Eq. ([!]) 

we get 

r = . (ii) 

(Cln(lO))^ 

Typical numbers of somatic hypermutations are N = 6 (Kiippers et al., 1993) correspond- 
ing to an affinity enhancement of 100 (£ = 2) (Janeway & Travers, 1997), while N = 8 
is considered to be a large number of mutations (Kiippers et al., 1993) corresponding to 
a high factor of affinity enhancement of 2000 (£ = 3.3) (Janeway & Travers, 1997). A 
huge affinity enhancement by a factor of 30000 (£ = 4.5) is reached with N = 9 mutations 
(Wedemayer et al., 1997). As can be seen in Fig. FJ], these values for the affinity enhance- 
ment with growing number of somatic hypermutations are consistent with a Gaussian 
affinity function (rj = 2) with width of F w 2.8. It should be mentioned that a consistent 
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Figure 2: The width of the Gaussian affinity function as a function of the factor of affinity en- 
hancement for different numbers of somatic hypermutations. Vertical lines denote the positions 
of the used experimental values of affinity enhancement. One can observe that consistently for 
all three values the resulting width becomes approximately 2.8. 

value for the width of the affinity function can not be obtained for other integer values of 
rj. But one should be aware that the used affinity enhancements are very rough estimates 
and thus our result should be understood as a first guess - even if the compatibility with 
the experimental data is convincing. 
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2.5 Shape space dimension 

The shape space dimension D enters the model through the mutation term in Eq. (|3D 
and has to be chosen appropriately to our model. We already mentioned that somatic 
hypermutation is likely to occur randomly (Weigert et al., 1970; Radmacher et al., 1998), 
i.e. that the direction of mutation in the shape space is not governed by the position of 
the antigen. The term describing mutation in Eq. (|3]) is in analogy to a diffusion term in 
a continuous space as used in (Perelson & Wiegel, 1999). Let us assume for a moment 
that the development of the GC reaction is governed by a diffusion process and let the 
number of maximum point mutations be iV max . Then the diffusion process operates in 
a sphere around a seeder cell of radius iV max in the shape space. In order to match the 
position of the antigen with at least one centroblast, all positions in the sphere have to be 
reached by diffusion. This means that the total number of centroblasts N c present in the 
GC 3 days after immunization, i.e. when selection is triggered, should be of the order of 
magnitude of the number of points in the sphere. The number of points in the sphere grows 
exponentially with the dimension of the shape space. Thus, for iV max = 8 and N c = 12000 
it follows that the shape space dimension should be 4. From the point of view of a pure 
diffusion process this is an upper bound for the shape space dimension. 

However, the above argument does not include the proliferation rate, the selection rate, 
and the recycling probability. The effective number of centroblasts available to reach the 
antigen may become larger than the 12000 cells assumed above through proliferation and 
recycling. As we want to describe a decreasing population for large times without further 
assumptions, the effective rate of production of additional centroblasts during the GC 
reaction cannot exceed some small value, since the centroblast population would explode 
otherwise for large times. So it seems unlikely that the shape space dimension appropriate 
for our model is substantially larger than 4. We will start with D = 4 and check our results 
for the dimensions 5 and 6. We want to remind that we are dealing with a phenotype space 
in this model and that a genotype space probably requires a substantially higher dimension. 



2.6 Long term behavior 

The parameter g was introduced in the model as the rate of differentiation from centrob- 
lasts in the dark zone into antibody presenting centrocytes ready to undergo a selection 
process. As we do not consider the centrocyte dynamics during the multi-step selection 
process (Lindhout et al., 1997), this parameter effectively describes the total duration 
of the selection process including the differentiation mentioned above, the selection at the 
FDCs and with T-helper cells, and the recycling to centroblasts. Only those cells that have 
successfully finished this program are contributing to the centroblast population described 
inEq.(|). 

Thus, being an effective measure of the speed of the selection process, we are able 
to deduce an upper bound for g. Solely the inhibition of apoptosis during the selection 
process at FDCs and with T-helper cells takes at least 1-2 hours (Lindhout et al., 1995; 
van Eijk & de Groot, 1999). This gives us an upper bound for the rate g of 
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Nevertheless, this upper bound is not sufficient to fix g. As g governs the whole reaction 
speed it will play a crucial role for the final state 21 days after immunization. As we know 
that at the end of the GC reaction only a few cells remain in the environment of the FDCs 
(Liu et al., 1991; Kelsoe, 1996), the parameter g will be tuned in order to get this result. 

Mathematically, this requirement has the form of a long term condition. Let us look 
for an equilibrium condition at the point 0o representing the antibody type of optimal 
complementarity to a presented antigen: 

dB , , . . 

— (0 = 0o,t) = . (13) 

Then, the selection term in Eq. (|3|) becomes simply (1 — q)ga B(<f) ,t) and for a non- 
vanishing centroblast distribution at 0o we get 

s = p-2pm-g + —^ + (l-q)ga = , (14) 
where we have introduced e as a measure for the vicinity to the equilibrium condition and 

B((po,t) 

/3(t) is a measure for the sharpness of the centroblast distribution in the shape space 
around the antigen position <p . (3 becomes constant when the selection process starting 
at day 3 after immunization has reached a phase in which the centroblasts at the position 
of the antigen dominate. This is already the case 8 days after immunization according to 
experiments (Jacob et al., 1993) and is confirmed in our model (see Fig.|]). 

Thus, at large times j3 can already be considered as constant in the equilibrium condi- 
tion Eq. (|TJD and the same applies to e. In this stage of GC development the centroblast 
population at 4> has a pure exponential form 

B(M = e £t . (16) 

The sign of e determines whether the centroblasts population explodes (e > 0) or dies out 
(e < 0). Because of the finite life times of GCs, e should adopt values slightly smaller than 
zero. 



2.7 Recycling probability 

The recycling probability 1 — q determines the fate of the positively selected centrocytes 
in the light zone. Apoptosis was inhibited for these cells and they may either become 
plasma or memory cells or return to the dark zone to reenter the fast proliferation phase. 
This recycling hypothesis (Kepler & Perelson, 1993) has been intensively discussed and 
the main position is that random somatic hypermutation - and it is likely that somatic 
hypermutation occurs randomly (Weigert et al., 1970; Radmacher et al., 1998), i.e. that 
the direction of mutation in the shape space does not depend on the position of the antigen 
- does not lead to a sustained optimization of affinity in a one-pass GC reaction (Oprea 
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et al., 2000). To reach a specific position in the shape space a multi-step optimization 
and verification process is required to avoid an arbitrary aimless walk through this high 
dimensional space. 

To propose a new perspective onto this question we consider the most convincing ex- 
perimental evidence for the existence of such a recycling process (Han et al., 1995b). A 
GC reaction was initiated in mice with one antigen. 9 days after immunization, i.e. in a 
stage of the reaction in which the good cells already dominate (Jacob et al., 1993), a second 
different, but related antigen was injected. This change of the reaction conditions in the 
course of the GC reaction leads to some very interesting observations: The rate of apop- 
tosis of centrocytes is enhanced by a factor of about !/«5 compared to the non-disturbed 
GC reaction. This can be explained by the fact that the antigen distribution in the shape 
space was changed such that most of the cells present after 9 days may fit to the first 
antigen but not to the second one. The probability of a positive selection is diminished 
and the inhibition of apoptosis does not occur to the same extent as before. On the other 
hand the total cell population in the GC vanishes with a u ~ 8-fold higher speed compared 
to the unperturbed GC reaction. This is very difficult to explain without recycling: As we 
know that the centrocytes (only very small numbers of centroblasts) undergo apoptosis if 
not rescued in time, the change of the apoptosis rate is of relevance only for the centrocyte 
population. The faster reduction of the whole GC cell population can not be explained by 
the enhanced apoptosis rate, if there exists no recoil effect of the centrocyte apoptosis to 
the centroblast population. This is an indirect confirmation of the recycling hypothesis. 

We want to consider this argument in a quantitative way by translating this experiment 
into the language of our model. We know from our consideration in Section PTB] that at day 9 
after immunization the optimization process is already completed and that the centroblast 
population at the antigen position = O in the shape space behaves according to Eq. ([TC|). 
But due to the modification of the antigen distribution at day 9 after immunization the GC 
reaction returns back to its dynamical selection phase. Thus, until day 9 the centroblast 
population Bx(<p, t) evolves with respect to the antigen distribution Ai(0*) as in Eq. (Q) and 
after day 9 two GC developments have to be compared: One continues its usual reaction 
as before. The second B 2 (4>,t) starts with the distribution Bi(<ft,t = 9d) and evolves with 
respect to the new antigen distribution 



where A0* denotes the shift of the additional second antigen with respect to the first one. 
Pi ;2 take into consideration the (possibly) different concentrations of both antigens. 

Now we are able to calculate the apoptosis rate in both scenarios for the cells of some 
fixed but arbitrary type 0. The number of centrocytes available for selection at every 
moment is gB(cj)). The apoptosis rate is given by this number reduced by the selected cells 



A 2 (4>*) 



Pi8(<t>* - 0o) + M(0* - (00 + A0*)) 



(17) 




(18) 



where i denotes the two scenarios. For the original GC reaction this becomes 




(19) 
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while in the GC with the new antigen distribution we find 



V 2 (<j>) = gB{(f>) (l-a o a(0-0 o )pi (! + «(</>))) 



(20) 



where we have shifted 0* in the second term of Eq. fllTD and introduced the ratio 

p 2 a(cj) - 0o - A0*) 



a(<p) 



pia(0 - O ) 



(21) 



valid for each point separately. The factor v of the apoptosis enhancement found in the 
experiment above is then given by the ratio of both apoptosis rates 



V2 _ 1 - a o a(0 - o )pi (1 + a(<j))) 
Vi 1 - a a((f) - 0o ) 



(22) 



On the other hand, this enhancement of apoptosis leads to a faster population reduction 
by a factor lu, which in our model is given by the ratio of the changes of the cell distributions 
over the shape space at the time of the presentation of the new antigen. According to 
Eq. ([!]) and using the new antigen distribution A 2 this ratio is given by 



dB 2 
dt 



dBi 
~dt 



-\ -1 



,t = 9d) 



p-2pm-g + if /?(0) + (1 - q)ga a{(t) - O ) px (1 + a(0)) 
p - 2pm - g + 22 /?(0) + (1 - q)ga a((f) - O ) 



(23) 



where the cell distributions -8(0, t = 9d) cancel as they are equal in both scenarios for 
t = 9d. 

Taking the results for v (Eq. (|22|)) and for uj (Eq. (|23|)) at the shape space point of the 
primary antigen = O the affinity function becomes equal to one. (3 is already a constant 
for t — 9 d such that we may eliminate pi(l + a(0o)) from both equations to get 



v 



p-2pm-g + ^ /3(0 O ) + (1 - q)ga 



uj-1 

or for the recycling probability 

1-q 



;i - q)g{ao - i] 



p - 2pm -g + ^(3((f) 



(24) 



a + (1 — a 



u-l 

UJ-1 



(25) 



This relation allows us to determine the parameter q and in this way to give a quantitative 
prediction for the recycling probability for selected centrocytes to reenter the dark zone 
and to continue proliferation. The recycling probability has to fulfill Eq. (|25|) in order to 
be consistent with the experiment of Han et al., (1995b). We want to emphasize, that the 



condition Eq. ( p5[) is independent of the concentrations of the two antigen types used in 
the experiment. 
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2.8 Output dynamics 



The production of plasma and memory cells optimized with respect to the presented antigen 
in the course of the GC reaction is governed by Eq. (||). According to our discussion of the 
different phases of the GC reaction in Sec. |272| , the starting time of output cell production 
is not necessarily correlated with the starting time of the mutation and selection process 
about 3 days after immunization. There may occur a time delay of the output production in 
relation to the selection process (see Fig. [[]). This possibility is encountered by comparing 
the dynamics of output cell production with experiment (Han et al., 1995a). Here, the 
number of output cells of high affinity to the presented antigen <fi is compared at day 6 
and day 12 after immunization and their relation is found to be approximately 

_ fw d dtO{<p ) _ 6 m) 
° edtO(<f> ) ~ 

As output cells are already present at day 6 after immunization, the starting point of output 
cell production has to fulfill Oh < At2 < 72 h (At 2 denotes the time window between the 
start of mutation and selection and the start of output cell production, see Fig. |l|). It is 
intuitively clear that for large At 2 the speed of output cell production vo will be enhanced, 
for the cell population of the GC will already be peaked at the position of the antigen in 
the shape space. Thus, we are able to adjust At 2 to the correct output cell dynamics. 



3 Results and Predictions 



Parameter 


quantity = 


= value 


Proliferation rate 


vl 


ln(2) = 


= 1/(6/0 


Somatic hypermutation probability 




m = 


= 0.5 


Rate of differentiation of centroblasts to centrocytes 


g/H2)- 


= 0.352 


Dimension of shape space 




D = 


= 4 


Selection probability for optimal centrocytes 




a = 


= 0.95 


Output probability of selected centrocytes 




q = 


= 0.2 


Number of seeder cells at FDCs 


Z*B{<f>,t 


= o) = 


= 3 


Radius of B-cell activation around antigen 




Ro = 


= 8-10 


Time duration of phase I of GC reaction 




Mi = 


= 72 h 


Time duration of phase II of GC reaction 




At 2 = 


= 48 h 


Time duration of the whole GC reaction 


El 


i AU = 


= 504/i 



Table 2: Summary of all parameters of the model. They were determined by experimental 
data. For explanations and references see the last section. The parameters g, q, and At 2 are 
determined in the course of the solution of Eq. (|) with respect to experimental data discussed 
before. 



Starting from a randomly chosen antigen epitope and - following the oligoclonal char- 
acter of the B-cell population of a GC - three randomly chosen activated B-cells in a sphere 
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Figure 3: Time course of the centroblast population integrated over the whole sphere in the 
shape space around the antigen (fi . In the first phase the initial B-cells proliferate. The large 
population reached at t = Oh is then reduced by the selection process. At t — 48 /i the 
production of plasma and memory cells starts, leading to an exponential decrease in the overall 
centroblast population. 



of radius Rq according to Eq. QUID, we let the GC reaction develop as described by Eq. (^) 
and Eq. (^j. The different phases (see Sec.^2|) of the GC reaction are respected. The 
parameters are determined from experimental data as summarized in Tab.|^. The param- 
eters g, 1 — q and At 2 , describing the differentiation rate of centroblasts to centrocytes, 
the recycling probability for the selected centrocytes and the time delay for the production 
of output cells, respectively, are iteratively adjusted in order to fit with the experiments 
described in the previous section. The time variable is running from to = — 72 h to accen- 
tuate the first phase of pure proliferation as a preGC-phase. The mutation and selection 
process is started at t — h and the production of plasma and memory cells begins at 
t = At2- The differential equation is solved numerically with a modified Euler method in 
a subspace, in which we assume Dirichlet boundary conditions with B = 0. 
The best results are obtained by an iteration procedure for 

q 0.355 . . 

, l-g = 0.8 , Atj = 48/i . (27) 



ln(2) h 

At t = h we find a large oligoclonal population of 12288 centroblasts as a result of the 
proliferation phase started with 3 initial cells (see Fig.|3|). These cells are of different type 
compared to the cells of optimal affinity to the presented antigen and in a typical distance 



17 



from them (varying between 3 and 8 point mutations to reach the optimal cell type). 

Optimization phase The selection process at the FDCs and with the T-helper cells 
starts and reduces the centroblast population with low affinity to the presented antigen. 
At the same time somatic hypermutation leads to a spreading of the centroblast distribution 
over the shape space, resulting in a small but non- vanishing amount of centroblasts at the 
shape-space position of the antigen (see Fig.[|). The large majority of B-cells does not 
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Figure 4: The centroblast population and the integrated plasma/memory cell production at the 
position of the antigen in the shape space. The centroblast population shows a growth due to 
proliferation and somatic hypermutation starting at t = Oh. The population growth is slowed 
down by the production of output cells starting at t = 48 h. Now the number of produced 
plasma and memory cells is increased steadily until the end of the GC reaction. Then, only 2 
cells (of optimal type) remain in the environment of the FDCs. 

survive the selection process (compare the summed cell population decrease in Fig. D), dies 
through apoptosis and is rapidly phagocytosed by macrophages present in the GC. 

Secondary proliferation As all positively selected centrocytes reenter the proliferation 
process in the optimization phase of the GC reaction, the number of centroblasts at the 
antigen position grows. Note that e > follows for the equilibrium condition Eq. ( |T4"D be- 
cause of q = 0, according to Eq. (f|). As a consequence, the overall centroblast population 
would restart to grow if this growth was not inhibited by the production of plasma and 
memory cells after 48 h. Especially, this secondary proliferation process is found for the 
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centroblasts encoding antibodies of high affinity to the antigen (see Fig. The centrob- 
lasts at the antigen position are practically all recycled (beside the reduction due to a < 1) 
such that no reduction process exists any more. This gives rise to a new perspective on 
the GC reaction: We need just as many initially proliferated cells, that a spreading of the 
distribution through somatic hypermutation leads to at least one or two cells matching the 
antigen position in the shape space. Then a considerable amount of cells of this optimized 
type is reached within the time delay At 2 of the production of output cells, opening the 
possibility of a full proliferation process for these cells. 

Recycling probability When the output production is turned on at t = 48 h the dy- 
namics of the germinal center is characterized by e < in Eq. (|T4"D, such that an exponential 
decrease of the whole centroblast population is initiated - including the population of the 
optimized cells (see Fig.|| and Fig.(|). Nevertheless, this decrease is slow enough to allow 
the production of a considerable amount of plasma and memory cells (see Fig. f| and Fig. |]). 
It is an interesting observation that a rather small number of centroblasts of the optimal 




Figure 5: A cut through the sphere in the shape space around the presented antigen at position 
2. The seeder cells are at a distance of at least 5 mutation steps. The number of produced 
plasma and memory cells grows considerably during the GC reaction. The rapidity of growth 
diminishes in course of time because of the decreasing number of optimized centroblasts at the 
antigen position. 

type with respect to the antigen is sufficient to produce this large number of plasma and 
memory cells - even with a small output probability of q = 0.2. On the contrary, a larger 
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output probability may decrease the integrated number of plasma and memory cells, as 
less cells reenter the proliferation process. In conclusion, we claim that a large recycling 
probability of 80% of the positively selected centrocytes is not only necessary to achieve 
an optimization by multi-step mutations but also to achieve a high number of resulting 
plasma and memory cells. 

The end of the GC reaction Finally, at day 21 after immunization, only a few cells 
with maximum affinity to the antigen remain in the environment of the FDCs, as required 
from experiment (Liu et al., 1991; Kelsoe, 1996). So the dynamics of our model allow 
for a vanishing of the germinal center population without any reduction of the antigen 
concentration or other additional requirements. The length of the germinal center reaction 
is determined by the interplay of the centroblast (into centrocyte) differentiation rate g 
and the delay of output production At 2 , which is determined by Eq. (|2~6"D , so that the GC 
life time is basically controlled by g. So from the perspective of our model a diminish- 
ing cell population may result solely from the interplay of proliferation, differentiation of 
centroblasts to centrocytes, and output production. 

Dependence of initially activated seeder cells The total number of remaining cells 
depends on the distance of the initial B-cells from the antigen <p . For numbers of mutation 
steps (necessary to reach the antigen) between N = 3 and N = 8, the number of remaining 
B-blasts varies from 43 to 0. For a typical number of mutations of iV = 5 about 10 B- 
blasts remain in the environment of the FDCs at day 21 after immunization. In other 
words, the GC life time depends on the initial distribution of the seeder cells with respect 
to the presented antigen. It would be interesting to check this statement experimentally 
by observing GCs with a variable number of mutations occurring during the optimization 
process. One could imagine for real GCs a self-regulating process which depends on the 
already positively selected centrocytes and thus prolongs the GC reactions with large 
distances between seeder cells and antigen, and shortens correspondingly the ones with 
small distances. However, there is experimental evidence that the intensity of the GC 
reaction depends on the quality of the initially activated B-cells (Agarwal et al., 1998). 
Also in quasi-monoclonal mice (Cascalho et al., 1996) it was found that the volume of GCs 
4 days after immunization depends on the average affinity of the reservoir B-cells to the 
antigen (de Vinuesa et al., 2000, Fig. 6). A quantitative comparison of these data with 
our model results would be interesting. 

3.1 Result stability 

The parameter set in Eq. ( p?[) is determined such that the experimental conditions Eq. fl25|) 
and Eq. (^) are fulfilled for typical initial conditions. Therefore, the prediction of the 
recycling probability to be 80% is a statement which - in the framework of our model - is 
justified to the same degree as the experiments are exact. Nevertheless, the sensitivity of 
the predicted high recycling probability to such modifications is weak. 

It should be mentioned that in addition there is a certain freedom of choice concerning 
the number of remaining cells at day 21 after immunization. For example, the parameter 
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set with a still higher recycling probability 

g 0.52 



ln(2) 



h 



q 



0.9 



Ato = 55 h 



(2i 



leads to a result of comparable quality as under the conditions of Eq. fl2~7|) - it is clear 
that a higher recycling probability slows down the production of output cells such that 
simultaneously the time delay of the production start has to be larger - with the difference 
that the number of remaining cells is reduced with respect to the previous result such that 
no cells remain at the end of typical GC reactions. For this reason and because g is on its 
upper bound Eq. flT2]), this parameter set was not considered as reasonable. 
If one requires a smaller recycling probability the parameter set 
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ln(2) 



0.286 
h 



l-g = 0.7 



At 9 = 42 h 



(29) 



leads to correct values for v and Eq. (^) is also fulfilled. The reduced recycling probability 
accelerates the output cell production such that it has to be decelerated by a reduction of 
the time delay At 2 . Even if the differentiation rate stays in a reasonable range, the number 
of remaining cells in the GC at day 21 after immunization is about ten-fold with respect 
to Eq. (^7]). More generally one may consider the number of remaining cells for typical GC 
reactions with an average number of mutations necessary to reach the antigen. 



Recycling probability 1 — q 


0.5 


0.6 


0.7 


0.8 


0.9 


£^B(0,* = 21d) 


1476 


515 


134 


10 






The mentioned experimental data are all reproduced for each value of 1 — q, solely the 
number of remaining cells is left open. To ensure that this number is in accordance with 
the experimental statement that a a few proliferating B-blasts remain at the FDCs at 
this stage of GC development, we are led to a recycling probability of 80%. We want to 
emphasize that this result is determined very clearly, as for the other recycling probabilities 
the remaining number of cells in the GC final state differs by at least an order of magnitude. 

This variation of parameters shows the range of stability of our prediction: 
First, the window of possible time delays At 2 , which are in accordance with Eq. (|26"D is very 
narrow. Already for At 2 > 55 h the condition Eq. ( |I~2"D will be violated. On the other hand 
for At 2 < 38 h, despite the immense number of remaining cells in the final GC state, the 
condition Eq. ( p6|) can not be fulfilled anymore. In consequence, it is not possible that the 
production of plasma and memory cells already starts with the establishment of light and 
dark zone, i.e. that it starts with the existence of a fully developed GC. The emphasized 
experimental bounds require necessarily a time delay for the production of output cells of 
at least 42 h. This means that the production process should be initiated by a separate 
signal and is not present in a working GC automatically. 

Secondly, the window of possible recycling probabilities is strongly bound by the experi- 
ments (Han et al., 1995a, 1995b) and by the experimentally observed number of remaining 
cells at day 21 after immunization. This leads to the conclusion that the recycling proba- 
bility should be at least 70% and not larger than 85%. Even ignoring the required number 
of cells in the final state, Eq. (|26|) does not allow a recycling probability smaller than 60%. 
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We would like to emphasize that this result differs extremely from the value for the recy- 
cling probability of 0.15 assumed in (Kesmir & de Boer, 1999). Nevertheless, we do not 
regard this discrepancy as a contradiction because in the model of Kesmir & de Boer not 
all parameters were fixed by experimental data. 

Shape space dimension The analysis presented above was performed with a shape 
space dimension of 4, (compare Sec. |2.5| ). It is an important remark that our results remain 
essentially unaffected for D = 5 or 6. Using for instance D — 5, only the parameter 
g, controlling the selection speed, has to be reduced slightly by 7% in order to stay in 
accordance with the experimental bounds. The total number of remaining cells at day 21 
after immunization does not change but its distribution is spread out on the shape space. 
For D — 6, the parameter g is again reduced by 6%, while all other parameters remain 
unchanged to stay in accordance with the experimental data. Again, the cell distribution 
at day 21 after immunization becomes less centered at the antigen position. We conclude 
that, for the values tested, the dependence on the shape space dimension is weak enough 
such that our results are not affected by it. Even, the slower selection speed for higher 
dimensions of less than 10% per dimension remains in the range of accuracy of our results. 

Note, that the dependence on the shape space dimension is restricted to an intermediate 
stage (phase II) of the GC reaction. In the first phase of pure proliferation D does not enter 
the dynamics at all. In the late phase of selection, when the distribution B(<p) becomes 
approximately spherically symmetric around the point representing the antigen, D cancels 
exactly with the factor 1/D in Eq. (£J) as there are D equal terms contributing to the sum 
such that D factorizes. A dependence of our results on the shape space dimension may 
only occur during the early phase of selection, in which the distribution B((p) is clearly 
not symmetric around the position of the antigen. Since the distribution of centroblasts 
develops into a spherically symmetric distribution around 0o already before day 8 after 
immunization (see Fig. ^), we deduce that a dependence on the shape space dimension may 
occur between day 3 and day 8 of the GC reaction only. 

Mutations with large changes of affinity We have assumed so far that mutations 
are represented in the shape space by next-neighbor jumps. To check the robustness of the 
result for the recycling probability against the possibility of mutations that lead to large 
distance jumps in the shape space, we consider a worst case argument. Let us suppose that 
10 % of the mutations lead to a large jump in the shape space. This is translated into the 
model by reducing the number of next-neighbor-mutations by a factor of 0.9. The result 
is a lower efficiency of the mutation process as the probability for a random mutation to 
improve affinity with respect to the antigen is very small. Therefore, the final number of 
cells is reduced by a factor of 10, and the centroblast to centrocyte differentiation rate has 
to be adjusted (from 0.244 to 0.255) in order to fulfill condition Eq. (^). A similar result 
as before is obtained by an adjustment of the recycling probability from 0.8 to 0.78. As this 
is a worst case argument, we consider our result as robust against large jump mutations. 
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Figure 6: The time course of (3((p ,t) (see Eq. (|15|)) at the shape space position of the antigen 
for the parameter set Eq.fl27|). (3 becomes constant shortly after t = 96 h, i.e. at day 7 after 
immunization. 

3.2 Optimization speed 

From experiment it is expected that the cells with high affinity to the antigen dominate 
already at day 8 after immunization (Jacob et al., 1993). This is verified in the model by 
checking if the asymptotic regime of the GC reaction is reached at day 8 in the sense that 
f3 (defined in Eq. ([15])) becomes constant. As can be seen in Fig. |6], (3 becomes constant in 
the course of day 7 after immunization for the parameter set Eq. Q2"7|). This behavior of 
(3 is not altered in a great range of parameter variation. Besides the reproduction of the 
experimental evidence for the optimization process to be accomplished after 8 days, this 
result is an a-posteriori check for the derivation of the recycling probability, which is based 
on a constant value of /3(<f>o) at day 9 after immunization (see Eq. ([24])). 

3.3 Optimization quality 

It may be interesting to verify if larger or smaller recycling probabilities lead to better 
optimization of the antibody types. The statement that recycling may be necessary to 
achieve a considerable affinity enhancement in a multi-step mutation process is widely 
discussed (Oprea et al., 2000) and enforced by the present work. We would like to point out 
that the relative number of good outcoming plasma or memory cells increases with larger 
recycling probabilities (see Fig. [7]). In the shape space this corresponds to a sharper peak 
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Figure 7: Comparison of the total plasma and memory cell distribution at day 21 after immu- 
nization on a cut along one of the coordinates in the shape space with the antigen at coordinate 
position 2. The parameters are uniquely determined by the experimental data with the excep- 
tion of the remaining number of cells at the end of the GC reaction, which is undetermined 
to allow different recycling probabilities 1 — q. The distribution of the output cells becomes 
sharper for larger 1 — q. 

of the distribution dtO((f),t). A too small recycling probability leads to a spreading 
of the outcoming plasma and memory cells in the shape space, i.e. to a weaker specificity 
of the GC output. 

On the other hand, we have already seen that a larger recycling probability lowers the 
absolute number of produced plasma and memory cells such that we are confronted with 
two competitive tendencies. Large recycling probabilities lead to specific but weak GC 
reactions, while small probabilities lead to unspecific but intense GC reactions. One may 
agree, that the value of 80% calculated here and supported by experimental data, is a good 
compromise between these two extremes. 

3.4 Start of somatic hypermutation 

There is experimental evidence that somatic hypermutation does not occur during the first 
phase of proliferation of centroblasts in the environment of the FDCs (Han et al., 1995b; 
Jacob et al., 1993; McHeyzer- Williams et al., 1993; Pascual et al., 1994b). However, we 
checked if our model allow somatic hypermutation to start during the proliferation phase. 
The primary reason for the retardation of the output production is to give time to the GC 
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to develop sufficient good centroblasts before the GC population is weakened through an 
additional output. Therefore, one may in principle think of a mutation starting long before 
the selection process. In this scenario the output production would start simultaneously 
with selection and one is led to a two phase process only: A first one with proliferation and 
somatic hypermutation and a second one with additional selection and output production. 
The centroblast-types would already be spread out in the shape space when selection 
and output start. There is no parameter set in accordance with all experimental data. 
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Figure 8: The centroblast population and the integrated plasma/memory cell production at 
the position of the antigen in the shape space. Except for the time phases we used the same 
parameters as in Eq. (0). The centroblast population shows a smooth enhancement due to 
somatic hypermutation starting at t = — 48/i. The population growth is stopped by the 
production of output cells starting at t — Oh. The slope of the summed output is substantially 
smaller than expected from experiment. 21 days after immunization only 1 cell remains. 

Especially, it is impossible to get the correct output dynamics as described in Eq. (^6|). 
Typically we get Vq = 2.9 being strongly different from the required value and a much 
smoother result (compare Fig.|8| to Fig. ^). In addition the distribution of plasma and 
memory cells is not as well pronounced at the antigen position such that the GC reaction 
gives rise to output cells with weaker specificity. Thus, we deduce from our model that 
the dynamic of plasma and memory cell production in real GC reactions does not allow 
somatic hypermutation to occur considerably before the selection process is started. 
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4 Conclusions 



We developed a new model for the GC reaction using its main functional elements. The 
GC reaction is described by the evolution of a centroblast distribution on an affinity shape 
space with respect to an initially fixed antigen distribution. On this shape space an affinity 
function, modeling the complementarity of antibody and antigen, was defined and its 
functional behavior as well as its width were deviated from affinity enhancements known 
for real GC reactions. 

In the model the reaction is decomposed into three phases (see Fig. |TJ) . In the pro- 
liferation phase a few seeder cells multiply. Our model excludes somatic hypermutation 
to occur already at this stage of the GC development. An optimization phase follows, in 
which in addition to proliferation, mutation of the antibody type and selection take place. 
This leads to a competition of spreading and peaking of the centroblast distribution in 
the shape space. In this phase, all positively selected cells reenter the proliferation phase 
in the dark zone, i.e. no plasma or memory cells are produced. This occurs solely in the 
output production phase. In this third phase all elements of the GC reaction are active. 

The dynamical evolution of the cell population in a GC is described by a set of coupled 
linear differential equation. All parameters are determined with experimental data and 
especially using (Han et al., 1995b). The model shows the typical behavior of GCs without 
any further adjustments. Due to the oligoclonal character of GCs a large amount of 
identical centroblasts is produced in the proliferation phase. These diffuse over the shape 
space in the second phase, in which only those cells are retained that exhibit a high 
affinity to the presented antigen. As the good cells reenter the proliferation process, this 
leads to a dominance of good cells within 7 to 8 days after immunization in accordance 
with experiments. The already optimized cell distribution then gives rise to plasma and 
memory cells leaving the FDC environment in the output production phase, leading to 
large numbers of output cells of optimized type due to the long duration of this third 
phase. The GC reaction is rather weak 21 days after immunization and only a small 
number of cells remains in the environment of the FDCs. This shows that the end of a GC 
reaction is not necessarily coupled with an antigen consumption or additional signals, but 
may occur simply due to the (unchanged) dynamical development of the GC. 

Our main result is that the experimental data and especially the evaluation of (Han 
et al., 1995b) lead to the following view of the GC reaction. The optimization phase lasts 
for not less than 42 h and not longer than 55 h. During this time all selected centrocytes 
are recycled and reenter the proliferation phase. Without such a non-output phase the 
number of optimal cells with respect to the presented antigen is not large enough to allow 
a considerable output rate during the remaining life time of the GC. We looked for an 
experimental evidence for this delayed production of output cells. Indeed, comparing 
two experiments (Jacob et al., 1993; Pascual et al., 1994a) we found that plasma and 
memory cells were first observed more than two days later than mutated cells during a GC 
reaction. This confirms experimentally the existence of the non-output phase that we found 
in the model. As output cells are not produced from the beginning of the GC reaction, 
the output should be triggered by a special signal, related to the affinity of the B-cells 
and the corresponding antigen. But even in the output production phase, the recycling 
probability turns out to lie between 70% and 85% and thus is substantially larger than 
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values expected until now (Kesmir & de Boer, 1999). We want to emphasize that such 
large recycling probabilities turned out to be crucial for a GC outcome of high specificity. 
Larger recycling probabilities lead to specific GC reactions of low intensity, while smaller 
ones to less specific but more intense GC reactions. The value of about 80% that we found 
represents a good compromise between these competitive tendencies. On the other hand it 
is not surprising that these recycled cells have not been found experimentally: The absolute 
number of recycled cells is very small compared to the total number of cells present in the 
GC, so that a measured signal for recycled cells is suppressed by more than two orders of 
magnitudes. A better signal should be observed in later stages of the GC reaction when 
the cells of high affinity to the antigen are already dominating. 

Furthermore, we conclude that the whole selection process including differentiation of 
centroblasts to antibody presenting centrocytes, the multi-step selection process itself and 
finally the recycling take around 3 or 4 h. This result is compatible with the experimental 
observation that the inhibition of apoptosis during the selection may take about 1 or 2 h. 

It is interesting, that the produced cells possess a certain broadness in the shape space, 
i.e. not only the centrocytes of optimal affinity to the antigen differentiate to plasma and 
memory cells. This may be of relevance for the resistance of an immune system against 
a second immunization with a mutated antigen. A considerable number of memory cells 
with high affinity to the mutated antigen may still be present for this reason. 

The results of our model suggest new experiments. A systematic and quantitative 
analysis of the GC reaction in dependence of the initial conditions, i.e. the number of 
necessary point mutations to reach the cell type of optimal affinity to the antigen, would 
lead to new insights into the interdependence of the essential elements of the GC. This may 
yield hints about self-regulating processes. Furthermore, it would be interesting to analyze 
quantitatively the total number of cells in the course of the GC reaction. Finally, it may 
be interesting to compare our model with its continuous counterpart. The translation of 
Eq. (|]) into continuous space leads to a differential equation of Schrodinger type with a 
Gaussian potential. 
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